Quantum Markovian activated surface diffusion of interacting adsorbates 
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A quantum Markovian activated atom-surface diffusion model with interacting adsorbates is pro- 
posed for the intermediate scattering function, which is shown to be complex-valued and factorizable 
into a classical-like and a quantum-mechanical factor. Applications to the diffusion of Na atoms on 
fiat (weakly corrugated) and corrugated-Cu(OOl) surfaces at different coverages and surface tem- 
peratures are analyzed. Quantum effects are relevant to diffusion at low surface temperatures and 
coverages even for relatively heavy particles, such as Na atoms, where transport by tunneling is 
absent. 
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In 1954, van Hove^ introduced the space-time corre- 
lation function G (a generalization of the well-known 
pair-distribution function from the theory of liquids) as 
a tool to study the scattering of probe particles off quan- 
tum systems consisting of ensembles of interacting par- 
ticles. Within the Born approximation in scattering the- 
ory, the nature of the scattered particles as well as the 
details of the interaction potential are largely irrelevant. 
Hence, following Lovesey^ the scattering processes with 
interacting particles essentially reduce to a typical prob- 
lem of statistical mechanics. The linear response func- 
tion of the interacting particles, also known as dynamic 
structure factor or scattering law, is then related to the 
spontaneous- fluctuation spectrum of such particles (mea- 
sured from G) by the fluctuation-dissipation theorem and 
can be expressed in terms of particle density-density cor- 
relation functions. 2 In general, G is a complex-valued 
function, what can be considered as a signature of the 
quantum nature of the problem. The imaginary part 
of G is important at small values of time (of the order 
of H(3, with f3 = 1/fceT), extending its range of influ- 
ence by decreasing the temperature. This dynamical 
regime takes place when the thermal de Broglie wave- 
length Xb = h/y / 2mkBT (m is the adsorbate mass) is 
of the order of or greater than the typical interparticle 
distances. The quantum system is then assumed to obey 
the fundamental condition of stationarity and the scat- 
tering problem satisfies the detailed balance principle, 
displaying the recoil effect. Here, we study the quantum 
observable effects of Na-atom diffusion on flat (weakly 
corrugated) and corrugated-Cu(OOl) surfaces probed by 
He atoms at different values of the Na coverage (9) and 
the surface temperature. For simplicity, in our model 
only coupling to substrate phonons (phonon friction) and 
not to low-lying electron-hole pair excitations (electronic 
friction) is considered. Nevertheless, as an extension of 
this stochastic model, the electronic friction could be ac- 
counted for by simply adding it to the total friction co- 
efficient. Moreover, diffusion by tunneling will not be 
considered in this work. 

The observable magnitude in this type of scattering 



experiments is the so-called differential reflection coef- 
ficient, which gives the probability that the He atoms 
reach a certain solid angle Q with an energy exchange 
Huj = Ef — Ei and wave vector transfer parallel to the sur- 
face AK = K/ — Ki after probing the quantum system, 
which consists of an ensemble of interacting Na atoms on 
the surface. This magnitude reads as 
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where n<j is the concentration of Na atoms on the surface, 
T is the atomic form factor, which depends on the He- 
Na interaction potential, and 5(AK, u>) is the dynamic 
structure factor or scattering law, which gives the line 
shape and provides a complete information about the 
dynamics and structure of the adsorbates. Experimental 
information about long-distance correlations is obtained 
from the dynamic structure factor for small values of AK, 
while information on long-time correlations is provided at 
small energy transfers, hu. 

The dynamic structure factor can be expressed as 



S(AK,uj) = J—jj e *( AKR - wt ) G(R, t) dtdR, 
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e~ lut I(AK,t)dt, 
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where R denotes the adatom position on the surface and 

/(AK, t) = ( e -iAK.R(0) e iAK.R(^ (3) 

In the Heisenberg representation, the trajectories R(i) 
are replaced by linear position operators R(t) and, there- 
fore, the brackets in Eq. ((3]) will denote the average of 
the expectation value of the operators enclosed. The 
so-called intermediate scattering function, /(AK, t), is 
the space Fourier transform of the van Hove G-function. 
Within the context of 3 He spin-echo experiments, this 
function is also known as polarization^- its real and imag- 
inary parts being observable magnitudes [for instance, in 
Cs diffusion on a corrugated Cu(001) surface^]. An ex- 
act, direct calculation of /(AK, t) is difficult to carry 
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out due to the noncommutativity of the particle position 
operators at different times. Nonetheless, this calcula- 
tion can be alternatively performed after postulating^ 
that the adatom position operators R(t) obey a standard 
Langevin equation (i.e., the associated stochastic dynam- 
ics is Markovian) and the product of the two exponen- 
tial operators in Eq. §3§ can then be evaluated according 
to a special case of the Baker-Hausdorff theorem (disen- 
tangling theorem), namely e A e B = e A+B e^- A ' B ^ 2 , which 
only holds when the corresponding commutator is a c- 
number. As shown below, this expression allows us to 
interpret /(AK, t) as the product of a classical-like and 
a quantum-mechanical intermediate scattering function 
associated with the exponentials of A + B and the com- 
mutator [A, B], respectively. Expressing now the thermal 
average implicit in Eq. ([3]) and the corresponding posi- 
tion operators in terms of the system Hamiltonian, and 
replacing t by t + ih(3, the so-called detailed balance con- 
dition reads as£ 



hods, its formal solution being 



S(—AK,—lj) 
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S(AK,u). 



(4) 



Recently, it has been shown^ that results obtained 
from a standard Langevin equation with two different 
non-correlated noise sources compare fairly well with the 
experimental data available. These noise sources are a 
Gaussian white noise accounting for the surface friction 
and a white shot noise replacing the pairwise interaction 
potential which simulates the adsorbate-adsorbate colli- 
sions. In this context, the double Markovian assump- 
tion holds because the substrate excitation time scale 
is much shorter than the characteristic times associated 
with the adatom motion (the maximum frequency of the 
substrate excitation is around 20-30 meV and the char- 
acteristic vibrational frequency of the adatom is about 
4-6 meV). Moreover, the time involved in a collision pro- 
cess is shorter than the typical time between two con- 
secutive collisions. Thus, memory effects can be ne- 
glected. Within this framework, called the interacting 
single adsorbate (ISA) approximation, the total friction 
n thus consists of the substrate friction 7 and the col- 
lisional friction A (i.e., n = 7 + A). The collisional 
friction can be related — to the surface coverage 9 as 
A = (6p0/a 2 )y/k B T/m, where a is the length of an as- 
sumed surface lattice with square unit cells and p is the 
effective radius of an adparticle. Thus, the stochastic 
single-particle trajectories R(i) running on the surface 
are assumed to obey the Langevin equation 



R(t) = -TjR(t) + F(R(t)) + 0N(t), 



(5) 



where F(R(t)) is the adiabatic force per mass unit de- 
rived from the periodic surface interaction potential, 
SN(t)=8NG(t)+5Ns(t) is the two-dimensional fluctua- 
tion of the total noise acting on the adparticle (the G 
and S subscripts stand for Gaussian and shot, respec- 
tively). In the Heisenberg representation, Eq. ((SJ still 
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- / §(nt-rf) F(R(t')) + <5N(i') dt',(6) 
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where R and P Q are the adparticle position and momen- 
tum operators at t = 0, respectively, and <j>(x) = 1 — e~ x . 
From Eqs. ([3]) and (|6|), and considering the disentangling 
theorem, J(AK, t) can be expressed as 



I(AK,t)~/,(AK,t)/c(AK,t), 



(7) 



which is a product of a quantum intermediate scattering 
function, J g (AK,t), and a classical-like one, / C (AK, t). 
The quantum contribution, governed by the commuta- 
tor, is the same for any type of surface regardless of 
its relative corrugation. For weakly corrugated surfaces, 
which can be assumed as flat, the commutator between 
the position and momentum operators is a c- number, 
and the relation (0 is exact. However, in general, the 
presence of the adiabatic force introduces an additional 
commutator, [R ,F(R(t))] = (ih)dF(R(t))/dP , where 
the dependence of the adiabatic force on the initial state 
(Ro, Pq) is through R(t). Assuming a Markovian regime 
(fast memory loss on the initial conditions), the previ- 
ous commutator is going to be negligible. Thus, in both 



cases, I q will reads as 



I q {AK,t) ~ exp 



i/iAK 2 
2nm 



${r/t) 



exp 



iE r <Z>(r]t) 



n 



(8) 

where E r = h 2 AK 2 /2m is the adsorbate recoil energy. As 
expected, the argument of I q becomes less important as 
the adparticle mass and the total friction increase (and, 
therefore, as the coverage also increases). Furthermore, 
for h — 0, we recover the standard classical function: 
I = I c . In ©, the time-dependence arises from Q(i]t). 
Thus, at short times (< h(3), we find $(ryt) « nt and the 
argument of I q becomes independent of the total fric- 
tion, increasing linearly with time. On the contrary, in 
the long-time limit, this argument approaches a constant 
phase. 

For flat or weakly corrugated surfaces, I c reads as 



ic(AK,f) 



-x a fat-*(i?t)] 



(9) 



where the so-called shape parameter \ is given by \ 2 = 
AK 2 (vq)/77 2 and, therefore, the total intermediate scat- 
tering function (jHJ) can be expressed as 



/(AK,*) 



D -X 2 lvt+aHvt)] 



(10) 



where a = 1 + ihrj/iksT if the thermal square veloc- 
ity is (vq) = ksT/m. Equation (Ti"Q|) is the generaliza- 
tion of the intermediate scattering function for the quan- 
tum motion of interacting adsorbates on a flat surface. 
The dependence of this function on AK 2 through the 
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FIG. 1: (Color online) Quantum intermediate scattering func- 
tion [Eq. (| 10 p ] for Na diffusion on a flat surface at 50 K and 
100 K: (a) real part and (b) imaginary part. Two coverages 
are considered: 9 = 0.028 (black solid line) and 9 = 0.18 
(red/dark grey dashed line). 



shape parameter is the same as in the classical theory^ 
Note that the total intermediate scattering function is- 
sued from Eqs. (fT|)- (fTU)) is exact for Gaussian quantum 
processes and no information about the velocity auto- 
correlation function is needed. However, classically the 
intermediate scattering function is usually obtained from 
Doob's theorem, which states that the velocity autocor- 
relation function for a Gaussian, Markovian stationary 
process decays exponentially with timeAi Furthermore, 
two regimes are clearly distinguishable: free-diffusion or 
ballistic and diffusive. The former is dominant at very 
low times, ijt -C 1, while the latter rules the dynamics 
at very long times, ijt 3> 1. In the diffusive regime, the 
mean square displacement, (R 2 (i)), is also linear with 
time, the slope giving the diffusion coefficient accord- 
ing to Einstein's law, D = ksT /mrj (which insures that 
the adparticle velocity distribution becomes Maxwellian 



asymptotically). Therefore, since tunneling is absent, the 
quantum diffusion coefficient follows Einstein's law as in 
the classical case. 

In Fig. [TJ the real and imaginary parts of Eq. ([7]) for 
Na diffusion on a flat surface are plotted at two different 
surface temperatures, 50 and 100 K, and two coverages, 
0.028 and 0.18. As can be clearly seen, the real part of 
I{t) decreases faster with temperature and slower with 
coverage. On the other hand, the imaginary part of I(t) 
displays maxima between 4-6% of the corresponding real 
part, depending on the temperature. It starts linearly 
with time and, after passing through a maximum, de- 
cays smoothly to zero. The corresponding quasielastic 
line shapes (around the zero energy exchange) will then 
display narrowing with the coverage and broadening with 
the surface temperature. This behavior could be experi- 
mentally confirmed for those systems where the diffusion 
barrier is smaller than the thermal energy ksT. For light 
particles, the imaginary part is expected to be much more 
important keeping the same shape. 

Let us now consider the case of nonzero corrugation. 
For Na atoms, the pairwise interaction potential is repul- 
sive and the mean interparticle distance should be greater 
than Xb most of time. Thus, the classical-like part of 
Eq. ([7]) could be replaced, at a first approximation, by 
the classical counterpart. Obviously, if diffusion is me- 
diated by tunneling, this procedure is no longer valid. 
The error comes from small times, but since diffusion 
is a long-time process, the influence on the quasielastic 
peak (wave-vector dependence) and the quantum diffu- 
sion constant (Einstein's law) will be very small for mas- 
sive particles. In Figs. [2] and 02 plots for Na diffusion on 
Cu(001) at the same surface temperatures and coverages 
as in Fig. [JJ are shown. The numerical values correspond- 
ing to I c have been obtained from classical Langevin 
simulations in a nonseparable adsorbate-substrate inter- 
action potential*^ The global time behavior of the real 
and imaginary parts of the intermediate scattering func- 
tion [see Figs. [U[a) and[3]Ja), respectively] display impor- 
tant differences for high and low temperatures. At 50 K, 
the thermal energy is not enough to surmount the diffu- 
sion barrier of the Cu(001) surface and adsorbates will 
remain for relatively long times inside potential wells. 
On the contrary, at 100 K, the intermediate scattering 
function decreases asymptotically to zero according to 
a more or less exponential function^ as happens in a 
regime characterized by surface diffusion. The maxima 
displayed by the imaginary part are, again, around the 
same percentage as before [see Fig. E^b)]. In principle, 
these imaginary parts should also be observable at least 
at 100 K in spin-echo experiments. With coverage, this 
time behavior means that the quasielastic peaks, given 
by the scattering law, will undergo broadening, as ob- 
served experimentally. 13 Quantum jump mechanisms can 
also be extracted from that peak in a similar way to the 
classical procedure^ The rapid oscillations displayed by 
the intermediate scattering function at short times [see 
Figs, dlb) and EJb)] indicate bound motions inside the 
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FIG. 2: (Color online) Real part (a) of the quantum inter- 
mediate scattering function [Eq. 0] for Na diffusion on a 
corrugated-Cu(OOl) surface at 50 K and 100 K. An enlarge- 
ment at short time-scales is shown in part (b). Two cover- 
ages are considered: 9 = 0.028 (black solid line) and 9 — 0.18 
(red/dark grey dashed line). 
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FIG. 3: (Color online) Imaginary part (a) of the quantum 
intermediate scattering function [Eq. 0] for Na diffusion on 
a corrugated-Cu(OOl) surface at 50 K and 100 K. An enlarge- 
ment at short time-scales is shown in part (b). Two cover- 
ages are considered: 9 = 0.028 (black solid line) and 9 — 0.18 
(red/dark grey dashed line). 



potential wells, which become less pronounced as temper- 
ature increases. These oscillations are associated with the 
lowest frequency mode or frustrated translational mode. 

In Fig. |H we show the effects of the quantum correc- 
tion in the diffusion process studied here at two different 
surface temperatures for a coverage of 0.028. To com- 
pare, both the classical intermediate scattering function 
and the real part of its quantum analog are displayed 
in the figure. As seen, although the Na atom is a rela- 
tively massive particle, at low temperatures the plateau 
is lower for the quantum case. This implies an initially 
(relatively) faster decay arising from the strong influence 
of the quantum behavior at short time scales. It is there- 
fore the real part of the intermediate scattering function 
what one should compare to the experiment rather than 



I Cl as is usually done. Obviously, this effect will be less 
pronounced at high coverages. 

It is remarkable that, within the Markovian approach 
presented here, the quantum intermediate scattering 
function, I q , is independent of the relative corrugation 
of the surface and, at short times, also independent of 
the friction. At low surface temperatures, the I q factor 
will be responsible for a higher contribution of the imag- 
inary part of J, given by Eq. ([7|) , modifying substantially 
the response in the diffusion process. Despite we have 
termed I c the classical-like intermediate scattering func- 
tion, strictly speaking it is not a classical magnitude, 
because it involves position operators. Nevertheless, for 
relatively heavy particles and at very long times (diffu- 
sion time scales), operators can be replaced by variables, 
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FIG. 4: (Color online) Classical intermediate scattering func- 
tion for Na diffusion on a corrugated-Cu(OOl) surface at 50 K 
and 100 K (black solid lines) and the real part of its quantum- 
mechanical analog (J7J (red/dark grey dashed lines). An en- 
largement at short time-scales is shown in part (b). The sur- 
face coverage considered here is 0.028. 

since As is very small. As far as we know, an exact quan- 
tum calculation for a corrugated surface is not possible 
and some approximations have to be invoked, e.g., the 
damped harmonic oscillator, which has been applied by 
some of usi^ within the same context. Of course, other 
different, alternative theoretical approaches can also be 
found in literature (see, for instance, Ref. [l5l ) within the 
single adsorbate approximation. The theoretical formal- 
ism that we propose here should also be very useful to 
avoid extrapolations at zero surface temperature when 



trying to extract information about the frustrated trans- 
lational mode. Diffusion experiments at low tempera- 
tures are very difficult to perform (or even unaffordable) . 
However, the type of theoretical calculations needed in 
this formalism is easy to carry out and they would pro- 
vide a simple manner to go to very low temperatures with 
quite reliable results, thus allowing to extract confident 
values of magnitudes such as friction coefficients and os- 
cillation frequencies. By decreasing the surface temper- 
ature, quantum effects are extended at higher values of 
time. Going from 100 K to 50 K, the time where the 
quantum dynamics is important increases from 0.07 ps 
to 0.15 ps. The standard propagation time for diffusion 
is greater than 400 ps as can be seen in Figs. [2] and [3] In 
our opinion, the range of applicability of this quantum 
theory should be around or below 10 K and with cover- 
ages up to 20%. Clearly, tunneling-mediated diffusion, 
where the isotopic effect and the so-called crossover tem- 
perature have been first observed by Gomer et al^ is not 
accounted for by our model. This is a very important as- 
pect which deserves further consideration, in particular, 
regarding the new observations of Zhu et al*^ concerning 
the coverage dependence of tunneling diffusion and the 
works by Ho and coworkers^ and Sundell and Wahn- 
strom. 19 

Finally, it is well known that the broadening is due to 
the increase of coverage. As we have shown in Ref. Hol . 
better agreement with the experimental results is found 
with the ISA model (classical theory) than with previ- 
ous calculations where the repulsive lateral interaction 
is taken into account in Langevin molecular dynamics 
simulations. This lead us to conclude that a stochastic 
description of the broadening should be good enough to 
describe it since the statistical limit (central limit theo- 
rem) in the number of collisions should be reached in the 
very long time propagation describing the diffusion pro- 
cess. This conclusion is still valid in this quantum theory 
since the adparticle is massive and the small quantum ef- 
fects are mainly observed in the short time limit. What 
this quantum theory has showed is the complex charac- 
ter of the intermediate scattering function whose real and 
imaginary parts are observable^ even for massive adpar- 
ticles (Cs atoms). 
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